Exact solutions to the spin-2 Gross-Pitaevskii equations 
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We present several exact solutions to the coupled nonlinear Gross-Pitaevskii equations which 
describe the motion of the one-dimensional spin-2 Bose-Einstein condensates. The nonlinear density- 
density interactions are decoupled by making use of the properties of Jacobian elliptical functions. 
The distinct time factors in each hyperhne state implies a "Lamor" procession in these solutions. 
Furthermore, exact time-evolving solutions to the time-dependent Gross-Pitaevskii equations are 
constructed through the spin-rotational symmetry of the Hamiltonian. The spin-polarizations and 
density distributions in the spin-space are analyzed. 
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I. INTRODUCTION 

The spinor Bose-Einstein condensates (BECs), which 
have been experimentally realized in optical potentials, 
exhibit a rich variety of magnetic effects. Due to the in- 
ternal degrees of freedom, they give rise to phenomena 
that are not present in the single component BECs^— , 
such as magnetic crystallization, spin textures, and frac- 
tional vortices, ect. In a spinor condensate, there exists 
an interplay between superfluidity and magnetism associ- 
ated to the spin-gauge symmetry in the Hamiltonian. A 
direct consequence is that the ferromagnetic BECs spon- 
taneously induce a supercurrent as the spin is locally 
rotated^ 5 -. 

The mean-field motion of the dilute spinor conden- 
sates is governed by the coupled Gross-Pitaevskii equa- 
tions (GPEs)£~— . There are many theoretical works that 
numerically solve GPEs or the corresponding stationary 
equations^—. Analytical solutions to the GPEs are gen- 
erally difficult because of the nonlinear density-density 
coupling between the atoms as well as the spin-spin cou- 
pling between the hyperhne states. Many efforts have 
been contributed to the one-dimensional (ID) soliton 
limit s 15 ' 16 , mainly for the F = 1 condensates. Never- 
theless, exact analytical solutions for the spinor BECs, 
especially for the F = 2 condensates, are absent in lit- 
erature. In a previous work, we have presented exact 
solutions to the F = 1 GPEs 17 . In this paper, we con- 
struct exact non-solitary solutions to the F = 2 GPEs 
in ID and give more insights into the internal structure 
of the states. The solutions are of complex form and 
are expressed in combinations of the Jacobian elliptical 
functions. 

The paper is organized as follows. In Sec. II we de- 
scribe the ID coupled nonlinear GPEs for F — 2 BECs. 
A set of generalized stationary equations are deduced. 
In Sec. Ill we present the special forms of analytical solu- 
tions in which only two of the components are nonzero. 
In Sec. IV, the associated time-evolving solutions to the 



special solutions are constructed by making using of the 
spin-rotational symmetry of the Hamiltonian. Section V 
includes a brief summary. 



II. EQUATIONS OF MOTION 

The spinor condensate formed by spin-2 atoms is de- 
scribed by a macroscopic wave function with five hyper- 
fine states ^ = (tp+2, ip+i, ipo, ip-i, i } -2) T - The mean- 
field Hamiltonian is expressed a o 18 i 19 

H = J dv { fit^V 2 + F(r)]V> m + 



m=— a 

|» a + ||P| a + ||^o| a }, 



(1) 



where F = ip^ n F m n 1 Pn with F (i = x, y, z) the 5x5 spin 
matrices. The coupling constants 5q, ci and c 2 are re- 
lated to scattering lengths Oq, a 2 and a 4 of the two col- 
liding bosons, with total angular momenta 0, 2 and 4, 
by c = 47r^ 2 (a 4 - a 2 )/7M, c x = 4nh 2 (3a 4 + 4a 2 )/7M, 
c 2 = 4-7r^ 2 (3a 4 - 10a 2 + 7a ) /7M. The total atom density 
isn = |^2| 2 + Ji/h| 2 + |?/'o| 2 + |V'-i| 2 + l'M 2 . The ampli- 
tude of the spin-singlet pair Aqq = (2?/> 2 -0_ 2 — 2tpiip-i + 
V>o)/v5- V(r) is the external potential. Hamiltonian (JTJ 
possesses the ^"(l)phase x ^'^(^)spin symmetry. The en- 
ergies are degenerate for an arbitrary state $ and its glob- 
ally spin- rotational states W — E/\f r , where U is the 5x5 
rotational matrix in the spin space which is expressed by 



the Euler angles as U(a,/3,'y) 



the ground state, the symmetry is spontaneously broken 
in several different ways, leading to a number of possible 
phase d 19 " 22 . 

We are concerned with the quasi-lD F = 2 BECs in 
a uniform external potential (V(r) = 0). The dynami- 
cal motion of the spinor wave functions are governed by 
idti/) m = SH/Sip^, which are explicitly written as the 
coupled nonlinear GPEs, 
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where F+ = Fl = 2(>|Vi + ^-1^-2) + V&(^o + 
r^-i) and F z = 2(|^ 2 | 2 - |^ 2 ) + |^| 2 - |V>_i| 2 . 
c o = co/2cij^, ci = c\j2a\ and c 2 = C2/2a\ are the 
reduced coupling constants with a± the transverse width 
of the quasi-lD system. 

Below we choose Ti = M = 1 as the units for conve- 
nience. Co, ci and c 2 are treated as free parameters. By 
substituting the wavefunction ty(x,t) with 



ipo(x,t) 



iptlx^e- 1 ^-^ 

^_ 1 {x)e~ i ^-^ t 
V V'-2(a;)e" 4(M+ ' Il)t / 



(3) 



we obtain the generalized stationary GPEs as, 



t2 



1-1% 

1 



c a n ± 2ciF z ]ip± 2 + ciF T ip ±1 



T2 



(p ± (ii)i/)±i = [--d x + c n±ciF z ]ip ±1 

+ Cl (^F T Vo + F ± ^ ±2 ) ~ ^ A ^T1 



1 /6 

-<9 2 + c n]^ + — ci(F + ^i + 



(4) 



It is notable that the solutions to the Eq.Q are not 
strictly "stationary" states since each component con- 
tains a distinct time-dependent phase factor, given /xi 
and H2 are not equal to zero. As we shall show, a 
"Lamor" procession in the spin space is associated to 
these states. The parameters /ii and /z 2 play the roles of 
linear Zeeman energies. Nevertheless, the density distri- 
bution of each hyperfine state are time-invariant. Hence 
we simply call these states as stationary states. 
The periodical boundary conditions 



iMl) = Vm(0), = C(0), 



(5) 



is adopted. In the following we consider two types of 
complex solutions to the stationary equations Q. 



III. SPECIAL SOLUTIONS 

In order to decouple the nonlinear spin-spin interac- 
tions in Eq.Q, we consider the special cases with only 
two of the hyperfine states are nonzero. The nonlinear 
density-density interactions are decoupled by making use 
of the unique properties of the Jacobian elliptical func- 
tions, as shown in the follows. The general forms of solu- 
tions in which all components are nonzero are obtained 
by the applying a rotation in the spin space, which will 
be addressed in the next section. 



A. Type A solution 

We first take the following ansatz as the solution to 
the nonlinear Eqs.fUJ), 



/ Mx) \ 

*Po(x) 
ip-x(x) 
W-2(a0 / 



i9(x) 






Dsn(kx, m) 




(6) 



where f(x) = yj A + Ben 2 (kx, m) and A, B and D are 
real constants, sn and cn are the Jacobian elliptical func- 
tions and k — AjK(m) with m the modulus (0 < m < 1). 
In the context we always take the number of periods j = 2 
as examples. One has 



sn 



D 2 



M-A 
B 



(7) 



By substituting the relations into the coupled GPEs(j4]) 
and making use of the relations between the Jacobian 
elliptical functions, one obtain the decoupled equations 



1 



^2^2 = ~2^' 2 ' +^ 2 IV'2| 2 V'2, 

M-i^-i = - 7^-1 +7-1 l 2 ^- 



(8) 



where the effective chemical potentials jl m and effective 
nonlinear coupling constants 7 TO are respectively defined 
as 



fa = n + H2 - c {A + D 2 ) + c n ^{B - D 2 ) 

-2 Cl (2A - D 2 ) + 2ci|(2B + D 2 ) 
72 = § (B - D 2 ) + ^-{2B + D 2 ), 



and 



fi-i=fi- fit - c (A + B) + 2cx{A + B) 

7-1 



- § I {D 2 -B) + §z{2B + D 2 ). 



(9) 



(10) 



The decoupled Eqs.flSJ can be self-consistently solved 
to yield, 
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FIG. 1: The density (a) and the phase (b) distributions for 
solution ([6]) with n = 2, m = 0.81. The physical parameters 
are Co = 8, a = 7, fx = 201.9003, m = 60, and fi 2 = 600. (c) 
The distribution of the spin-polarization. 



FIG. 2: Snapshots of the density distribution for ([6]) (attached 
by the time factor) in the spin space at t = 0, t = T/8, and 
t — T/3, respectively. The procession period is T = 2n/u). 



21.2 



5 



m k 

72 



.4 



2/2 2 - A: 2 (l - 2to 2 ) 



The phase in Eq.© 
0(x) 



372 

is calculated as 



P(0 



(12) 



(13) 



where ai = ±(2£l 2 A 2 - 2j 2 A 3 + k 2 (l - m 2 )AB)? is the 
integral constant. The periodical boundary conditions 
([5]) require that the amplitudes and phase satisfy 

/(l)=/(0), 0(1) - 0(0) = 2j7T x n, (14) 

respectively, where n is an integer. The periodical con- 
dition for the phase is satisfied by properly adjusting the 
modulus m of the Jacobian elliptical functions. In the 
calculations, we fix the modulus m and fix, [i 2 while ad- 
just /i to fulfill the periodical boundary conditions (fT4")) . 

Figure 1 display the distributions of the density, the 
phase and the spin-polarization of the solution ©. The 
relevant parameters are Cq = 8, C\ = 7, m = 0.81, 
fj, = 201.9003, /ii = 60, [ii = 600, and n = 2. Numeri- 
cally, we get jl 2 = 655.2336 and 72 = 24. It implies that 
the effective interactions in the m = — 1 hyperfine state 
should be repulsive while in the m = +2 hyperfine state 
should be repulsive. One notes that F x (x) = F y (x) = 
and A 00 (x) = 0. Obviously, the spin-polarization vector 
F(x) is not enough to exhibit the spin configurations for 
the F = 2 condensates. In order to reveal the "Lamor" 
precession of the spin, in Fig. 2 we display the density 
distribution in the spin space. At the nodes of the V'-i 



(x = 0, 0.5, and 1.0) it becomes fully spin-polarized with 
the spin-configurations exhibit the axial symmetry. At 
other spatial positions the spin-configuration exhibit the 
tctrahedral symmetry, which subsequently rotate an an- 
gle to the adjacent positions. As we have stressed, al- 
though state ([6]) is the solution to the stationary equa- 
tion [[%]). it contains a spin procession with frequency 
uj = \i% + Mi- Fig. 2 snapshots this procession at t = 0, 
t = T/8, and t = T/3, where period T = 2tt/uj. 



B. Type B solution 

We next take an alternative ansatz as the solution to 
the stationary Eqs.((3]), 

i6(x) 



( M*) \ / f(x)e^ \ 
ipi(x) 
ipo(x) 
i/t-i(x) 

V 1>-2(X) ) 



(15) 



\ Dcn(kx, m) ) 



where f(x) = y/A + Bsn 2 (kx,m) and A, B and D are 
real constants. One has 



l^2| 2 



(16) 



D 2 ' B 
By the same way, we obtain the decoupled equations 



as 



A-2^-2 



f//i 2 + 7-2^-2 |V 



(17) 
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periodical boundary conditions (|14p and the modulus m 
is accordingly determined. 




FIG. 3: The density (a) and the phase (b) distributions for 
solution (|15l) with n = 1, m = 0.4. The physical parameters 
are c = -100, a = -80, c 2 = -20, fx = 55.8083, /Lti = -40, 
and /J.2 = 20. (c) The distribution of the spin-polarization. 



FIG. 4: Snapshots of the density distribution for (|15|) in the 
spin space at t = 0, t — T/8, and t = T/3, respectively. The 
procession period is T = 2-k/uj. 



where the effective chemical potentials jl m and effective 
nonlinear coupling constants 7 m are respectively defined 
as 

M2 = M + Mi - co04 + £> 2 ) + c |(S - D 2 ) 

-4 Cl (A - # 2 ) + 4 Cl |(B + £ 2 ) - *$j£-(A + B)(18) 



7 2 = C o(l-^)+4c 1 (l + ^)- 2 ff 



and 



£_ 2 = M + M2 - c (A + B) + 4ci(A + B) - if {A + B) 

B 

TP 



7-2=c (l-^)+4 Cl (l + ^)-^ 



(19) 

The decoupled Eqs. (jl"7|) are self-consistently solved to 
yield 



A-2 = ifc 2 (l-2m 2 ), 7_ 2 



L» 2 



(20) 



Figure 3 display the distributions of the density, the 
phase, and the spin-polarization of the solution (|15p . The 
physical parameters are chosen as Co = —100, c\ = —80, 
c 2 = -20, fj, = 55.8083, m = -40, and [i 2 = 20, n = 1 
to = 0.4. Numerically, we get //2 = 33.0909 and 72 = 
-208. Analogously, F x (x) = F y (x) = and A 00 (x) ^ 0. 
In order to exhibit the spin procession of the state, we 
show in Fig. 4 the density distribution in the spin space at 
t = 0, t = T/8, and t = T/3, respectively. The procession 
frequency is oj — (J,2 — Ml an d T = 2tt/uj. The spin is fully 
polarized as the nodes of hyperfine state ij)-2- At other 
positions, the density distribution has C^ z symmetry and 
is obviously distinct to that of solution © . 



B = 



i 2 k 2 
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A = 
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The phase in Eq.([T5| is 

6(x) = 



a 2 



(21) 



(22) 



where a 2 = ±(2fl 2 A 2 - 2j 2 A 3 + k 2 AB)^ is the integral 
constant. The amplitudes and phase should satisfy the 



IV. GENERAL SOLUTIONS 



The more general time-evolving solutions to the GPEs 
@ can be obtained by making use of the symmetry of 
the Hamiltonian. If \&(x, i) is a solution to the GPEs 
([2]), then 9'{x,t) = U*ff(x,t) is also a solution, where 
U is an arbitrary rotational transformation in the spin 

space which is explicitly written with the Euler angles 

"^23 
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(23) 



In the following we fix the three Euler angles as a — the time- evolving solutions to Eqs.([2]), 
7tt/13, (3 = tt/4 and 7 = tt/11. 



A. Solution associated to Type A 

By attaching the time factors to the solution ([6]) and 
then applying the rotational transformation, we obtain 



1// = e~ v * 



■j) 2 cos 4 | e -2i(a+7) e -V2' 
p 2 COS 2 I 
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(24) 



/ 



where 02 = ^/A + Bcn 2 (kx, m)e l6 ^ and 0_i = 
Dsn(kx,m). It should be noted that even at f = 0, the 
solution (|M|) is no longer a solution to the stationary 
Eq. Q . The temporal evolution of the density distribu- 
tion of each component is shown in Fig. 5. The rightmost 
column of Fig. 5 displays the snapshots of the spin rota- 
tion at x = 0.143. 



B. Solution associated to Type B 



Similarly, we obtain the general time-evolving solution 
associated to (fT5|) as, 
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(25) 



/ 



where </>i = \J A + Bcn 2 (kx, rnje lS ( x > and = 
Dsn(kx,m). Figure 6 show the temporal evolution of 
the density of solution (|25l) . The rightmost column is 
the snapshots of the spin rotation at x = 0.2. 



V. SUMMARY 



In summary, we have presented two classes of ana- 
lytical stationary solutions to the ID coupled nonlinear 
GPEs which govern the dynamics of the spin-2 conden- 
sates. Obviously, we can obtain a lot of other exact solu- 
tions with different combinations of the Jacobian ellipti- 



6 



cal functions. The exact time-evolving solutions are also 
constructed. The spin-polarization and the spin proces- 
sion are addressed. 
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FIG. 5: Temporal evolution of the density of solution (|24[l . 
The parameters are the same as in Fig.l. The rightmost col- 
umn displays the spin rotation at x = 0.143. 
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FIG. 6: Temporal evolution of the density of solution (|25|l . 

The parameters are the same as in Fig.3. The rightmost col- This work is supported by the funds from the Min- 
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